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ABSTRACT 

We present two Rossi Xray Timing Explorer (RXTE) PCA observations of the type 2 Seyfert 
galaxy NGC 4388 caught in an unusual low X-ray absorption state. The observations were trig- 
gered by a detection in the 1.5-3 keV band of the RXTE All-Sky Monitor. NGC 4388 was found 
at a somewhat high continuum level (f(2-10 keV)=8xlO _11 erg cm _2 s _1 ) and with a column density 
(Njj ^3xl0 22 cm -2 ), a factor ~10 lower than normal. The second PCA observation, four hours later, 
gave Njf <2xl0 21 cm -2 indicating, at the 3.1<r level, variability so rapid puts the absorber on a few 
100 Schwartzschild radii scale, similar to the Broad Emission Line Region, or smaller. This small 
scale creates difficulties for the parsec-scale obscuring torus paradigm of Unified Schemes for type 1 
and type 2 AGNs. 
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1. INTRODUCTION 

Optically NGC 4388 is a classical type 2 Seyfert galaxy 
(Huchra, Wyatt & Davis 1982) with permitted and for- 
bidden emission lines of the same width (Khachikian & 
Weedman 1974). There is abundant evidence that many, 
and perhaps all, type 2 Active Galactic Nuclei (AGNs) 
are normal type 1 AGNs with both the characteristic 
broad emission lines and the optical to X-ray continuum 
obscured by a flattened torus of absorbing gas and dust 
(e.g. Mulchaey et al. 1994). This is the basis of the Uni- 
fied Scheme for AGN (Antonucci 1993, Urry & Padovani 
1995). NGC 4388 has been detected in X-rays for over 20 
years (Table 1) and has always shown a column density, 
Njj-= 2-5 xlO 23 cm" 2 . The most common form of the 
Unified Scheme locates this absorption in a dusty torus 
at parsec distances from the central continuum (Krolik 
& Begelman 1988, Pier & Krolik 1992, 1993). However, 
Risaliti, Elvis & Nicastro (2002) found that 23/24 X- 
ray absorbed AGNs (10 22 cm -2 <N# <3xl0 23 cm -2 ) 
showed N# variability by a factor 2-3. The best stud- 
ied objects varied on the shortest accessible timescale of 
months, which is rather fast to be due to Keplerian mo- 
tion at parsec radii and so raises questions about the 
nature of the obscuring torus. 

Risaliti et al. suggested an alternative location in the 
cool outer parts of an accretion disk wind, echoing the 
model of Kartje, Konigl & Elitzur (1999) who predicted 
just such a torus. This location predicts much faster 
Nh variability, down to a timescale of days. In a simple 
model of Poisson variations in the number of obscuring 
clouds, N c , the amplitude of variability found by Risal- 
iti et al., implies N c ^5-10. In this case 0.1-1% of the 
time N c =0 and, in the Unified scheme, the central type 1 
nucleus would then be unveiled. 

The RXTE (Swank et al. 1998) All Sky Monitor 
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(ASM, Remillard et al. 1997) is just sensitive enough 
to detect such low energy 'unveiling events'. We thus 
began a Target of Opportunity (TOO) program with 
Rossi-XTE to obtain snapshot Proportional Counter Ar- 
ray (PCA, Swank et al. 1998) spectra of type 2 AGNs 
showing signs of a low energy detection in the RXTE 
ASM. Here we report the detection of a low N# 'unveil- 
ing event' in NGC 4388. 

2. OBSERVATIONS AND DATA REDUCTION 

We monitored NGC 4388 with the RXTE ASM to 
search for detections in the soft 1.5 - 3 keV X-ray 
channel ('a'). Normally NGC 4388 has a flux of 
~4xl0 -13 erg cm -2 s -1 in this band (Forman et al. 
1979), while a detection requires a flux some 250 times 
larger (~lxl0 _lo erg cm" 2 s _1 ). Simply removing the 
large absorbing N# (~2-5x 10 23 cm~ 2 ) would increase the 
observed flux to ~4xl0 _11 erg cm -2 s — 1 , a factor of 100, 
so that only a modest additional factor 2-3 increase in the 
emitted continuum would be needed to put NGC 4388 
over the threshold for ASM detection. By contrast, an 
increase in the emitted continuum flux by a factor >100 
would be unprecedented among the well studied type 1 
AGN, where factors of a few to ~10 variation are seen 
(Markowitz et al. 2003). Hence an ASM 'a' band detec- 
tion is a good indicator of a low Njj event. A trigger- 
ing event of this type occured on 2003 May 9 (Fig. 1), 
shortly after another one (which is visible on the left of 
Fig. 1). One day after the ASM trigger NGC 4388 was 
observed twice with the RXTE PCA, for 1.9 ksec and 
6.2 ksec, with a four hour gap between the two observa- 
tions (Fig. 1). 

We only consider data from PCU-2 for this analy- 
sis, as this is the best-calibrated proportional counter 
unit (PCU) in the RXTE/PCA 5 . Data reduction tools 
from LHEASOFT version 5.2 6 were used to screen and 
prepare the event files and spectra. Data were taken 
in "Standard 2" mode, which provides coverage of the 
full PCA bandpass (2-60 keV) every 16 seconds. Only 

5 URLs http://lhcawww.gsfc.nasa.gov/uscrs/kcith/bkgd_status/status.html 
http: / /hcasarc . gstc. nasa.gov / hstscrv / grodis/ msgUUUtifci . html 
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Fig. 1 — The RXTE All-Sky Monitor (ASM) light curve of NGC 
4388 which triggered the pointed observations. The 'a' (1.5-3 keV) 
band fluxes are shown as circles and the 'c' (5-12 keV) band fluxes 
are shown as triangles . The triggering data points are marked. 
The two pointed PCA fluxes in the same bands are also shown. 
The error bars on the PCA measurements are smaller than the 
points. 



data from the top Xe-filled gas layer in PCU-2 was 
used to make the source and background spectra, as 
this gas layer has the lowest background. The standard 
RXTE "good time interval" filtering was applied to the 
data. For the non-imaging PCA accurate background 
subtraction is crucial for faint sources like NGC 4388. 
The background spectra were made using the tool 
pcabackest using the latest 'faint source' background 
model (pca_bkgd_cmfaintl7_e5vv20031123.mdl, updated 
by RXTE in November 2003). pcabackest calculates 
the predicted (dominant) particle background every 16 s, 
and so tracks the variation of the background around the 
orbit, thus taking into account the different backgrounds 
in these two short observations. Redistribution matrix 
files (rmfs) and ancillary response files (arfs) were made 
and combined into a single instrumental response file us- 
ing pcarsp. 

We added 0.6% systematic errors to our spectra using 
grppha, as we find that in many instances acceptable 
fits to the Crab can be obtained with 0.6% systematic 
errors. However, Poisson errors of 5% - 10% are domi- 
nant. The lowest channels in each of the PCUs routinely 
reveal strong deviations likely due to calibration uncer- 
tainties; in addition, the calibration of the PCUs is more 
uncertain above approximately 25 keV, and the spec- 
tra of faint sources like NGC 4388 become background- 
dominated in this regime. In fitting the spectra, then, 
we ignored energy range below 3 keV (channels 1-4) and 
above 20.0 keV. The PCA X-ray spectra are shown in 
Fig. 2a, where they are compared with a Chandra ob- 
servation performed 14 months earlier. Fig. 2a shows 
that most of the variation is at low energies, <5 keV. In 
Fig. 2b we show the ratio between the two RXTE spec- 
tra, which indicates that the cut off in the spectra at low 
energies is due to a difference in Nfj. 

To fit the PCA spectra we used XSPEC v.11.3 7 . A 
model comprising: a power-law of slope T, absorbed 
by a zero redshift Njy, with a superimposed emission 
line near the 6.4 keV Fe-K line energy, was fit to the 
two PCA data sets and gave good \ 2 (30/34 degrees 
of freedom for PCA-1 and 25/34 degrees of freedom for 

6 URL http://heasarc.gsfc.nasa.gov/docs/software/lhcasoft/ 
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Fig. 2.— Upper panel: RXTE spectra of NGC 4388. For com- 
parison, we also plot the best fit model of the Chandra observation 
performed one year earlier, convolved with the response and effec- 
tive area of the PCA, and extrapolated from 10 keV to 20 keV 
for clarity. Lower panel: Ratio PCA-2/PCA-1, showing greater 
variability at low energies. No change in N^j would give a flat line. 



PCA-2. The low redshift of NGC 4388 (2535 km s" 1 
Huchra, Wyatt & Davis 1982) is indistinguishable from 
zero with the PCA. The Galactic N H (2.5xl0 20 cmT 2 , 
Murphy et al. 1996) is negligible for the PCA energy 
range. The total flux in the kiloparsec scale X-ray neb- 
ula around NGC 4388 (Iwasawa et al. 2003) is also neg- 
ligible (2.6x 10 _13 erg cm _2 s _1 ). The results are given 
in Table 1, together with fits to the same model, with 
the same minimum energy, for two unpublished observa- 
tions from XMM-Newton, and for six data sets from the 
literature. The measured N# in the two PCA observa- 
tions is ~ 5 x 10 22 cm -2 and < 0.9 x 10 21 cm" 2 , respec- 
tively (90% confidence). This rapid change between the 
two PCA observations in four hours is significant at the 
3.2(7. I.e. there is a 0.14% chance that the N# value is 
the same in the two PCA observations. 

The PCA Njf values are one and two orders of magni- 
tude respectively lower than in all previous X-ray obser- 
vations of this source. In particular the PCA-1 N# is 13ct 
smaller than that measured by XMM-Newton 5 months 
earlier (XMM-2, Table 1). [The change in N H in the 
5 months between XMM-1 and XMM-2 is significant at 
the llcr level.] The contour plot of Fig. 3 shows that 
the dramatic reduction in Njj is not due to a degener- 
acy between T and N# that is sometimes encountered in 
X-ray spectra. T is flat (r=0.8-0.9) in the PCA spectra 
compared both with most of the earlier observations and 

7 Arnaud K., & Dorman B., 

2003. XSPEC 11.3 User Guide, URL 

http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/manual.html 
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Fig. 3. — Photon index (r) vs. Column density (Nff) contour 
plots (68%, 90% and 99% confidence), derived here, for the two 
PCA observations of May 2003, and the two XMM-Newton obser- 
vations of July 2002 and December 2002. 



with unobscured type 1 AGNs (T —1.8, Nandra et al. 
1997, Reynolds 1997, Perola et al. 2002). The XMM-1 
observation (Table 1) also gave a flat T, yet is heavily 
obscured, so a flat spectrum is not a property correlated 
with a low Nff. The small (—10 arcsecond) beam size of 
XMM-Newton (Jansen et al. 2001) effectively rules out 
the possibility that the flat PCA spectral slope is due to 
another source lying within the — 1 sq.degree PCA field 
of view. The high 2-10 keV flux in the BeppoSAX-1 ob- 
servation demonstrates that a similarly high flux does 
not reduce the observed Nff. So we must be seeing bulk 
motion across our line of sight, as in Risaliti et al. (2002). 

The N# variation is robust against reasonable changes 
in the model: (1) Adding a reflection component 
(PEXRAV model, Magdziarz & Zdziarski 1995) does not 
alter T or Njj- Even leaving all the parameters free, 
the best fit is obtained with a covering factor (R=Q/2tt) 
R=0. The 90% confidence upper limits are R<0.8 for 
the first observation, and R< 0.2 for the second one. 
Fixing R to these upper limits we obtain column densi- 
ties N H i ~ 2 x 10 21 cm" 2 and N H2 - 8 x 10 22 cm" 2 . 
(2) Adding a soft component, of course, does change the 
fit dramatically as this component affects only the lowest 
channels. For black body emission at typical soft excess 
temperatures of kT —0.3-1 keV the required 0.1-2 keV 
flux is 20 mCrab, or a luminosity of 5x 10 42 erg s _1 . This 
luminosity causes no physical problems although a 10 4 s 
variation requires an optically thick source (Elvis et al. 
1991 ApJ 378, 537); a sphere of the implied area has a 
radius of only 10 10 cm, 0.3 light-seconds. However such 
a large low energy flux is not suggested by the lowest 
energy PCA bins that were omitted from the fit. 

The x 2 t es t prefers an emission line consistent with the 
6.4 keV Fe-K line, but broadened by cr=0.42t°^ 4 keV 
(90% confidence), (i.e. FWHM = 2.38<r = 1.0 keV, 
FWHM/E = 0.15c). A broad line would make an origin 
for the line in a parsec-scale torus unlikely, and would ar- 
gue for wind or accretion disk origin. The XMM-2 obser- 
vation however, using the higher spectral resolution CCD 
EPIC detectors, gives a lower gaussian a (73±20 eV), 
though the line may be complex. The PCA-2 measured 
equivalent width of the line, —500 eV, is stronger than 



normal for a type 1 Seyfert, but less strong than can be 
found in Compton thick AGN (Levenson et al. 2002). 
Again XMM-2 gives a more normal EW (-200 eV) for 
an AGN that is not Compton thick. The EW and a of 
the PCA and XMM-1 are consistent but both disagree 
with the XMM-2 values. 

3. DISCUSSION 

We have found a factor 100 decrease in the col- 
umn density toward the normally almost Compton thick 
(r es -0.1-0.3) type 2 AGN NGC 4388. This decrease 
certainly occured in less than the 0.4 years from the ear- 
lier XMM-Newton observation. The obscuring material 
in type 2 AGNs is thus in a highly dynamic state, and 
warrants intensive monitoring. A few strongly Comp- 
ton thick AGN have become almost Compton thin on a 
timescale of 2.5 - 5.5 years (Matt, Guainazzi & Maiolino 
2002), but not on shorter timescales, and still with resid- 
ual column densities of — 10 23 cm" 2 even in the low ab- 
sorption state, in 3 out of 4 cases. NGC 4388 is the 
only known case of an AGN in which a substantial X-ray 
opacity changes to an undetectable value (r es <0.001). 

Moreover, it is likely that the decrease in obscuring col- 
umn density coincided with the 2-day 'flare' seen in the 
RXTE ASM, so that the flare is primarily an 'unveiling 
event'. Without X-ray spectra through the rise of the 
ASM flux, however, we cannot be certain. Between the 
two PCA observations we saw evidence, at 3.2cr, for a 
change of N# in four hours. 

The probable short timescale (either —2 days or 
4 hours) of the column density variations has strong im- 
plications for the location of the obscuring matter. As- 
suming that the absorption is due to clouds in Keple- 
rian motion around the central source, and interpreting 
the 4 hour time lag between the two observations as the 
crossing time of a cloud implies a distance from the cen- 
ter R< 10 4 p 2 t|i?5, where p w is the density in units of 
10 10 cm -3 , Rs is the Schwartzschild radius, and t± is the 
timescale in units of four hours (Risaliti et al. 2002). If 
NGC 4388 has an Eddington ratio of 0.1 then the black 
hole mass is — 10 6 - 1O 7 M . Scaling from the K-band 
bulge magnitude gives a similar mass. This implies that 
the absorber is at a distance typical of the Broad Emis- 
sion Line 'clouds' (BELRs), or smaller, and is of similar 
density. Only if a high density (p > 10 12 cm" 3 ) is as- 
sumed can a parsec distant absorber produce changes 
on the observed timescale (Fig. 4). A similar conclusion 
applies to the Seyfert 1.5 Galaxy NGC 4151, from the 
detection of a A(Njj) — 2 x 10 23 in a time interval of 
~ 150 ksec during a BeppoSAX observation (Puccetti et 
al. 2003). 

This result is a challenge to the parsec scale usually 
attributed to the obscuring torus. For spherical isolated 
clouds, the high density implies small (r c — 10 10 cm) 
cloud sizes: (1) The expansion time for a 100 K (v sounc i — 
1 km s _1 ) cloud from 10 10 cm to 10 11 cm is only 10 6 s,, 
far shorter than the orbital time. But the large pressure 
needed to confine them (p/k — 10 14 K cm" 3 ) cannot 
be provided by self-gravity. (2) A hot confining medium 
would produce thermal emission much greater than the 
observed AGN bolometric luminosity (assuming a layer 
of 1 parsec thickness, L gas = 10 48 nf erg s" 1 , where 
T7 and 717 are the density and temperature in units of 
10 7 K and 10 7 cm" 3 , respectively, Rybicki & Lightman 
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Instrument N° T E b Fi , ev ^f k f<! l/ Obs. Date Rei h 

SL2-XRT 2.1ZH 1.9l'o 5 2.1 3.5 1985 Jul 29 T~ 

ASCA-1 3.15±i J 1.32to.?I 6.49t°;°| 150i^ 750±| 00 1.3 1.5 1994 Jul 04 2 

ASCA-2 3.34±o.9 l^Ttois? 6 - 47 ±ao7 < 200 700 - 250 °- 64 °- 7 1995 Jun 21 2 

BeppoSAX-1 3.80t°' j 1.58±°;22 ^-^-o.il < 230 233 ±35 5 2 5 3 - 5 1999 Jan 09 3 

BeppoSAX-2 4.80tJ' | lA7±Q °f 6.38lg;°g < 120 5251^2 0.94 1.4 2000 Jan 03 3 

Chandra 3 -50t°' j 1.8 j 6.36±o;°2 < 230 440tg° 0.36 0.6 2001 Jun 08 4 

Chandra-0 2.50±°' 2 1 - 25 -0.2S 6 - 36 ±a03 < 130 165 ±60 2 - 9 2 - 8 2002 Mar 05 ^ 

XMM-1 l-70±o 0.9lt°;37 6.4lt[j;°2 < 77 503tg° 0.77 0.7 2002 Jul 07 5 

XMM-2 2.6lj;o 1% l-46toog 6 - 44 to 02 73 ±20 204 ±26 2 -° 2 - 2 2002 Dec 12 5 

, + 0^5 n qq+0.11 a o.+O.'ll , oon r n , + 138 



RXTE/PCA-1 0.521^ 24 0.99Iq ■{{ 6.34+^J < 380 5031^ 7.2 4.0 2003 May 10» 5 

h0.07 a o 9 + 0:i0 qnn + iSO c;7n+114 
-0.03 °- ,: ' z -0.09 3au -190 °' u -100 



RXTE/PCA-2 < 0.09 0.86t° o3 6.321 390 -i90 570 ±ioo 6- 1 2 - 9 2003 Ma y 109 5 



TABLE 1 

a Absorbing column density, N h , in units of 10 23 cm -2 . b Peak energy of the iron Ka line, keV. c Width of the iron line, 
in eV. d Equivalent width of the iron line, in eV. e Observed 2-10 keV flux, in units of 10 -11 erg s _1 cm -2 , f Intrinsic 2-10 
keV luminosity, in units of 10 42 erg s _1 cm -2 (assuming a Virgo cluster location at 20 Mpc, Tammann et al. 1999). S XTE 
observations 4 hours apart, h References: 1: Hanson et al. 1990; 2: Forster, Leighly & Kay 1999; 3: Risaliti 2002; 4: 
Iwasawa et al. 2003; 5: This work. j Fixed. k 0-order spectrum from HETGS observation (Obsid=2983). 
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Fig. 4. — Distance (R/R s ) vs. density (p) plot for AGN X- 
ray absorbers. The "standard torus" of the unified model has a 
density n < 10 6 cm -3 . Variability of the absorber in 4 hours or 
less rules (solid diagonal line) out an absorber with such properties, 
and is more compatible with the standard parameters of the Broad 
Emission Line Clouds (BELR). 



1979). (3) The N# through a confining medium would 
be 10 25 5 nj x d pc cm -2 = 10r es . So all continuum vari- 
ations shorter than a few years would be smeared out by 
Thomson scattering, contrary to observations (Table. 1). 
(4) Individual clouds would cover only a small fraction, 
~10 -6 , of the X-ray emitting source, whose dimensions 
are > 10 13 cm for a black hole of mass > 10 7 M Q . If 
hundreds of clouds are needed to cover the X-ray source, 
no Nh variation by more than a few percent can be ob- 
served with significant probability. 



To produce large variations in N# compatible with < 
iV c >^5-10 needs clouds of a diameter within a factor 
of a few of the continuum source size. This implies a 
density of ~10 9 cm~ 3 . This is just compatible with a 
radius R ~ /ewl00i? s given a 4 hour variation (Fig. 4). 
If the absorbers are part of a wind crossing our line of 
sight (Elvis 2000) then sheet-like structures (Arav et al. 
1998) become plausible, and this allows larger radii. 

The great majority (>99%) of the obscuring gas in 
NGC 4388 occurs at small radii, and the most likely 
scenario is that Broad Emission Line clouds are drift- 
ing across our line of sight in NGC 4388 leading to large 
changes in N#. Similar behavior has been seen in some 
type 1 AGN, involving an N# ^10 23 cm~ 2 over ^100 d at 
the radius of the BELR (Lamer, Uttley & Miller 2003). 
Since we have no simultaneous optical spectra we do not 
know whether the bulk of the dust, which absorbs the 
optical and UV photons, lies on the same small scale. 
The closest dust could be to the continuum is the subli- 
mation radius (Barvainis 1987) which, in NGC 4388 at 
~10 42 erg s" 1 (2-10 keV), is -7xl0 16 cm (0.02 pc), or 
2xl0 4 r s (A/ 7 ) (fig.4). The dust:gas ratio in AGN is typi- 
cally a factor 10 below the Milky Way value (Maccacaro 
et al. 1982, Maiolino et al. 2001), so that a separate 
dusty absorber is allowed. Unveiling events in type 2 
AGN, such as the one reported here for NGC 4388, put 
strong constraints on Unified Models for AGN, and seem 
to point to a different view of the obscuring torus. 
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